Electron-Phonon or Hole Superconductivity in MgB27 
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The BCS electron-phonon mechanism and the unconventional 'hole mechanism' have been proposed 
as explanations for the high temperature superconductivity observed in MgB2- It is proposed that 
a critical test of which theory is correct is the dependence of Tc on hole doping: the hole mechanism 
predicts that Tc will drop rapidly to zero as holes are added, while the electron-phonon mechanism 
appears to predict increasing Tc for a substantial range of hole doping. Furthermore, the hole 
mechanism and electron-phonon mechanism differ qualitatively in their predictions of the effect on 
Tc of change in the B — B distances. We discuss predictions of the hole mechanism for a variety of 
observables as a function of doping, emphasizing the expected differences and similarities with the 
electron-phonon explanation. The hole mechanism predicts coherence length and penetration depth 
to increase and decrease monotonically with hole doping respectively. 



I. INTRODUCTION 

Superconductivity at 4QK in MgB2 was not predicted 
by theory. After its discovery [|l|, it has been pro- 
posed that this finding is expected within two fundamen- 
tally different theoretical frameworks: the BCS electron- 
phonon theory |^,^ and the theory of hole superconduc- 
tivity [Q. Both theories have claimed to be consistent 
with various experimental observations. The purpose of 
this paper is to expand on the predictions of the theory of 
hole superconductivity, and to make sharper the distinc- 
tion between it and the electron-phonon theory before 
critical experiments are performed that can differentiate 
between both theories. It is generally easier to differen- 
tiate between theories by comparing their predictions , 
rather than their postdictions, of experimental obser- 
vations; elaborate theoretical frameworks can often find 
consistent explanations even for the most unexpected ob- 
servations. 

The electronic structure of MgB2 is well established, as 
a variety of old ^ as well as new calculations are 

in essential agreement. Approximately 30% of the den- 
sity of states at the Fermi energy is due to planar boron 
Px,y states (cr bonds) that have little dispersion in the 
z direction, giving rise to nearly cylindrical hole Fermi 
surfaces of 2D character. The remaining 70% of the den- 
sity of states originates in boron pz states (tt bonds) that 
are strongly hybridized with the Mg s — p orbitals, have 
3D character and give rise to mostly electron-like Fermi 
surfaces. No d-electrons exist in either Mg ot B, so that 
magnetic and strong correlation mechanisms (generically 
called 'big tent' mechanisms [^) proposed for the high 
Tc cuprates, do not appear to be applicable. 

Within the electron-phonon framework, two different 
explanations have been proposed, hereafter referred to 
as EPl 1^ and EP2 [^. Both explanations emphasize 
the importance of strong bonding of the boron atoms in 



giving rise to strong electron-phonon coupling, as well 
as the light ionic mass giving rise to a large prefactor in 
the BCS-Eliashberg expression for the transition temper- 
ature. They appear to differ in the relative contribution 
of the boron states. Whereas EPl appears to suggest 
that contributions from all states are important, EP2 at- 
tributes superconductivity exclusively to the nearly full 
boron px,y states. It is argued that the observation of 
a boron isotope effect Q (isotope coefficient a = 0.29) 
strongly favors electron-phonon mechanisms j^J^,^ . 

In contrast, within the theory of hole superconductiv- 
ity (l^Jll[] the electron-phonon interaction is irrelevant, 
and instead superconductivity originates in undressing of 
hole carriers, driven by Coulomb interactions, in bands 
that are almost full. The superconducting condensation 
energy is kinetic, since paired carriers have lower effec- 
tive mass than unpaired ones, and electron-hole symme- 
try breaking is central to the physics. In MgB2, the fact 
that large parts of the Fermi surface are strongly hole- 
like, together with the fact that the boron planes where 
the holes propagate are highly negatively charged, are 
roposed to be the essential factors giving rise to high Tc 
The existence of an isotope effect is generically ex- 
pected within this theory also , although its magni- 
tude is difficult to calculate; a simple estimate yields Q| 
a much larger isotope effect than observed experimen- 
tally in MgB2 if the electron-phonon coupling suggested 
in is used. 

Both the electron-phonon theory and the hole theory 
predict that the superconducting state is s-wave, which 
appears to be supported by tunneling jl^-p^ as well as 
NMR ||l6| measurements, and both theories are consis- 
tent with the observation of an isotope effect Q . The the- 
ory of hole superconductivity requires Q that the con- 
ductivity in the normal state is hole-like, which is consis- 
tent with recently reported Hall effect measurements [0 ; 
it also requires pt] superconductivity to disappear when 
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the hole bands in MgB2 become full, which is consis- 
tent with reported experimental results on Mgi^xAlxB2 
psf . Both of these facts apparently are consistent with 
electron-phonon theory, as discussed in EP2 On the 
other hand, the hole theory predicts that Tc should in- 
crease under pressure if the dominant effect is reduction 
of the B — B distances [|j , while electron-phonon theory 
(EPl predicts that pressure generically should reduce 
Tc; exper iments show that hydrostatic pressure reduces 
Tc [ [l9[|2C[| , in apparent agreement with electron-phonon 
theory. We return to this point later in the paper. 

A key prediction of both theories is the behavior of Tc 
upon hole doping, for example in the compound (not yet 
fabricated to our knowledge) LixMgi^xB2- Here there 
is a clear opportunity for distinction between both the- 
oretical frameworks. EPl explicitly states that decreas- 
ing the Fermi level 'may provide an additional contribu- 
tion to A', which suggests that Tc should increase upon 
hole doping. EP2 does not explicitly address this crucial 
point, but emphasizes the 'substantial value of the Fermi 
level density of states', even though 'the hole density n/j 
is small'. Since the density of states increases with hole 
doping, both for the a as well as for the tt bands, for a 
very substantial range of hole doping, a reasonable in- 
ference within electron-phonon theory is that Tc should 
also increase with hole doping for a substantial range. 
Under the assumption that electron-phonon matrix el- 
ements, phonon frequencies and Coulomb pseudopoten- 
tial don't change substantially, electron-phonon theory 
would predict that Tc should increase with hole doping 
in a range of about 2eV below the Fermi level of MgB2 
(following the density of states increase) , corresponding 
to ^ 0.36 holes added per B atom, and stay high (above 
the Tc of MgB2) till about AeV below e_F, corresponding 
to ~ 0.61 holes added per B atom. Instead, the theory of 
hole superconductivity predicts that Tc will rapidly drop 
when holes are added, becoming small or zero before the 
number of added holes per B atom reaches only 0.12. 
Comparison between these predictions of both theories 
will be discussed in detail in the next sections. No mea- 
surements on the behavior of Tc with hole doping have 
yet been reported to our knowledge. 



width, U the on-site Coulomb repulsion, W is propor- 
tional to the nearest neighbor Coulomb repulsion, and K 
is proportional to the correlated hopping interaction At 
that drives superconductivity. The BCS pairing interac- 
tion is given by Vkk' = V{ek,ek'), with 



V{e,e') = U + 2^{e + e')+A^ee'. 



(1) 



The critical temperature is determined by the equation 
1 = 2Kh - Wh - UIo + {K^ - WU){Iol2 - Ii) (2) 

and the parameters and c that define the energy- 
dependent gap 



A(e) = A, 



Ti/2 



by the equations 

1 = K{h + do) - W{l2 + ch) 



c = K{l2 + ch)-U{Ii+ch) 



with 



D/2 



dJ_^/2'^'^ D/2^ 2E{e) 



E{e) = v/(7^7^j^TW 



(3) 



(4a) 



(4b) 



(5a) 



(5b) 



with / the Fermi function and /i the chemical potential. 
Finally, the hole density is determined by the equa- 
tion, 

2 f^^\ ^ l-2f{E{e)) 



D 



D/2 



2E{e) 



III. CHOICE OF PARAMETERS 



II. MODEL OF HOLE SUPERCONDUCTIVITY 

Within the model of hole superconductivity, as well 
as within EP2, the bands that drive superconductivity 
are the nearly two-dimensional boron pa bands. Calcu- 
lations for three-dimensional anisotropic band structures 
I pof have shown that a two-dimensional model reproduces 
the essential features, hence we will ignore the third di- 
mension here for the calculations with the model of hole 
superconductivity. We will also approximate the nearly 
constant density of states of the hole pa bands by a con- 
stant, which has a negligible effect. The model is then 
defined by 4 parameters: U, K, W and D. D is the band- 



The density of states at the Fermi level of MgB2 is 
estimated to be 0.75 states/eV, of which approximately 
0.25 states/eV is ascribed to the B pa states In a 

model with constant density of states and a single band, 
that would correspond to D = AeV . There are two pa 
bands that contribute to this density of states, approx- 
imately with 2/3 and 1/3 weight g (heavy and light 
hole bands respectively in the nomenclature of An and 
Pickett). According to An and Pickett, the equivalent 
two-dimensional flat bands have densities of states 0.18 
states/eV and 0.07 states/eV, which would correspond to 
bandwidths Di = h.&eV , D2 = \AeV respectively. The 
total number of holes in the B pa bands for Mgi?2 is esti- 
mated to be nt = 0.13/unit cell, of which approximately 
0.09 and 0.044 holes are in the heavy and light hole bands 
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respectively. The total rih per B atom in MgB2 is ap- 
proximately 0.067. 

We will discuss elsewhere the results of our theory in 
the presence of two hole bands, which we do not expect 
will be qualitatively different Here, we will use a sin- 
gle 'effective band' of bandwidth D = beV , Coulomb re- 
pulsion U — beV and nearest neighbor repulsion W — Q. 
As we will discuss in a later section, the model gives sim- 
ilar results for a wide range of parameters. For these 
parameters, we choose the value of K required to yield 
Tc AQK for uh = 0.067, which is K ^ 2.%7eV. As dis- 
cussed elsewhere |Q] , we believe this value of K is reason- 
able for MgB2, but emphasize that here K is a fitting pa- 
rameter as we have not obtained it from a first-principles 
calculation. 



IV. Tc VERSUS DOPING 

As discussed in the introduction, no calculations of Tc 
versus hole doping with the electron-phonon model (EPl 
or EP2) have yet been reported. It is possible that such 
calculations may indicate large changes in the phonon fre- 
quencies, electron-phonon matrix elements or Coulomb 
pseudopotential with doping. In the absence of other 
information however, we will assume that all these quan- 
tities are constant with doping, and calculate within 
the electron-phonon model by the modified McMillan for- 
mula 01 



< UJlog > -1.02 , '\+^'., 

1.2 



(7) 



with n* =0.1 and < ujiog >= TOOK |], and A propor- 
tional to the density of states in the B pa band. Be- 
cause the fractional contribution of the B pa bands to 
the total density of states is quite constant with dop- 
ing (1/3), this calculation should predict the results of 
both electron-phonon models EPl and EP2 under the as- 
sumption that electron-phonon matrix elements, phonon 
frequencies and /i*do not change with doping. 

Figure 1 shows the results of this calculation for the 
two models, as function of hole content in the pa bands 
per B atom, rih- Note that is approximately 1/3 of 
the total hole doping per unit cell, n^°*. The electron- 
phonon model predicts that Tc will continue to increase 
well beyond the point where the hole model predicts Tc 
will have vanished; the maximum Tc of 9iK occurs for 
n/i = 0.43 per B atom, or total hole doping of MgB2 of 
approximately ~ 1.1 holes per unit cell, correspond- 
ing to bringing the Fermi level down approximately 2.leV 
from its position in MgB2. In contrast, the maximum Tc 
in the hole model of 49if occurs for ^ 0.035 per B 
atom, corresponding to electron doping of MgB2 of ap- 
proximately 0.03 electrons per B atom, or 0.09 electrons 
per unit cell. 

As mentioned above, these results could be modi- 
fied in the electron-phonon model if there are substan- 



tial changes in some or all phonon frequencies, electron- 
phonon matrix elements, or Coulomb pseudopotential. 
In the hole model, some modification may be expected if 
the contribution of the two B pa bands is taken into ac- 
count separately. In particular, assuming it is the heavy 
hole band that dominates Tc, Tc would go to zero in an 
even narrower range of hole doping than indicated by Fig- 
ure 1. Despite these caveats, we believe the qualitative 
difference in the behavior predicted by the hole model 
and by the electron-phonon models depicted in Figure 1 
is robust. 



V. PRESSURE DEPENDENCE OF Tc 

In the hole model, a decrease in the B — B intraplane 
distances should strongly increase Tc- So far, only ob- 
servations of changes in Tc under hydrostatic pressure 
on polycrystalline samples have been reported |p^ , po| , 
that indicate that such pressure decreases Tc- We be- 
lieve that hydrostatic pressure is likely to affect much 
more strongly the lattice spacing in the c direction than 
the planar lattice spacings, due to the stiffness of the pa 
bonds. Furthermore, it is possible that substantial charge 
transfer occurs between different bands when pressure is 
applied. For example, in many high Tc cuprates the hole 
concentration in the planes is increased by approximately 
10% when IGPa hydrostatic pressure is applied ||2^. Ac- 
cording to Figure 1, a 10% increase in hole content from 
MgB2 leads to a decrease in Tc of iK; the reported ob- 
servation of a decrease in Tc by 1.6^ |l^ could hence 
be accounted for by such an increase in the B pa or- 
bitals hole content together with a small decrease in the 
B — B distances. Hence the observation is not necessarily 
inconsistent with our model, as also emphasized in ref. 

Here we consider the effect on Tc of a decrease in the 
B~B intra-plane distances. The intrinsic effect of chang- 
ing lattice spacing in the c direction should be much 
smaller within our model p^ . In the electron-phonon 
model, we assume the dominant effect will be to decrease 
the density of states. In the hole model, we assume the 
effect is to increase the bandwidth (i.e. decrease the den- 
sity of states) and increase the interaction parameter K , 
that depends on overlap matrix elements as the band- 
width does, by the same fraction. 

Figure 2 shows the changes expected under a 5% 
change in these parameters, achieved by either phys- 
ical or chemical pressure. In the electron-phonon 
model (again assuming no change in phonon frequen- 
cies, electron-phonon matrix elements and /i*), a small 
decrease in Tc results . In the hole model, a strong in- 
crease in the critical temperature for all hole dopings re- 
sults. By performing such experiments and monitoring 
the changes in lattice constants and in carrier concentra- 
tion (e.g. through Hall measurements) we hope it will be 
possible to decide which of the two qualitatively different 
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behaviors shown in Figure 2 takes place in this class of 
materials. 



VI. OTHER RESULTS FOR THE MODEL OF 
HOLE SUPERCONDUCTIVITY 

We next discuss other results for our model for a single 
band with constant density of states, in the range of pa- 
rameters that may be appropriate for this class of mate- 
rials. Figure 3 illustrates the effect of changing Coulomb 
interaction parameters in the model, always choosing the 
parameter K so as to yield the observed value Tc ^ iOK 
for Uh = 0.065. It can be seen that the behavior of 
Tc versus doping is quite insensitive to large variations 
in the Coulomb interactions. On increasing the near- 
est neighbor repulsion the range of hole dopings where 
Tc is nonzero increases somewhat, and on increasing the 
on-site Coulomb interaction that range decreases. If the 
Coulomb repulsion appropriate for AIgB2 is larger than 
5 eV, the maximum Tc obtained by electron doping could 
be larger than 50 K, as seen in Figure 3. 

Figure 4 illustrates the effect on Tc of changing the 
bandwidth, i.e. the density of states. It can be seen that 
the effect is again remarkably small, with a reduction in 
the bandwidth leading to a small decrease in the range 
of hole concentration where Tc is non-zero. 

We next calculate other observables for the parameters 
of Figure 4. Figure 5 shows the behavior of the coherence 
length ^o, defined as the average size of the pair wave 
function, with hole doping. The formulas to evaluate this 
quantit y fo r the model under consideration here are given 
in ref. [Q. The coherence length is found to be almost 
independent of the values of the Coulomb interactions, 
but it depends strongly on the bandwidth, as seen in the 
different curves in Figure 5: as the bandwidth decreases, 
the coherence length decreases. The coherence length 
in Figure 5 is given in units of lattice spacings in an 
effective square lattice; to transform to physical units 
for MgB2 (a = 3.14A), a lattice spacing Oe// = 2.22 A 
should be used. For hole concentration Uh = 0.065, the 
result obtained with D = 2.5 eV is close to the observed 
experimental value, = 45 A [ |25| , p6[ |. Instead, for D = b 
eV we obtain a coherence length of = 83 A, larger 
than seen experimentally. However our calculated value 
corresponds to the in-plane coherence length which is not 
necessarily the same as that measured in a polycrystalline 
sample. As seen in Figure 5, the coherence length is 
predicted to increase monotonically with hole doping. 

Similarly, Figure 6 shows the behavior of the in-plane 
London penetration depth Xl, assuming the clean limit: 

A,^ 4638(d(A))V2 1 

Ta[meV)\'-'^ 

with d the distance between boron planes, d = 3.52 A, 
and Ta the average in-plane kinetic energy per boron 
atom. The estimated value for MgB2 is Al = 1400 A 



[ |25[ , which is close to the value given in Figure 6 for 
D = 2.5 eV and Uh = 0.065, Xl = 1344 A. For D = 5 eV 
we obtain a smaller value than seen experimentally. The 
penetration depth is predicted to decrease monotonically 
with hole doping. In conjunction with the increasing co- 
herence length, this implies that the Ginzburg-Landau 
parameter k = Xl/^o will rapidly decrease with hole dop- 
ing, which could eventually lead to a cross-over from type 
II to type I behavior for high hole doping. However this 
is likely to be prevented by disorder, that would cause 
an increase in the penetration depth from its clean limit 
value. 

The gap versus hole concentration follows closely the 
behavior of the critical temperature. This is shown in 
Figure 7. In contrast to high Tc oxides, we do not find a 
substantial increase in the gap ratio in the underdoped 
regime, because the parameters here correspond to a 
weak coupling regime. In Figure 8 we show the hole con- 
centration dependence of the specific heat jump at T^, 
which agrees with the BCS weak coupling value 1.43 for 
high hole concentrations and becomes larger for low hole 
concentrations, particularly as the bandwidth becomes 
smaller. 

The temperature dependence of various quantities ob- 
tained from our model also follows closely the BCS weak 
coupling behavior. As an example we show results for the 
gap ratio and the specific heat for one parameter set in 
Figure 9. Experimental results for specific heat of A/gi?2 
show a clear specific heat jump at the transition with 
value close to the expected BCS value p7|| . 

Finally Figure 10 shows tunneling characteristics for 
one set of parameters and hole doping appropriate to 
MgB2. Again the behavior resembles the weak coupling 
BCS results for an s-wave gap, except for the existence of 
asymmetry. As emphasized elsewhere pO| , an asymmme- 
try of universal sign occurs for this model, with a larger 
peak for a negatively biased sample. The case of Figure 
10 corresponds to the smallest bandwidth considered; for 
larger bandwidth the magnitude of tunneling asymmetry 
decreases. 



VII. CONCLUSIONS 

The growing number of experimental results on Af gi?2 
suggests that superconductivity in these materials is 
more akin to conventional superconductivity than it is 
to high temperature superconductivity in the cuprates. 
Thus it is natural that a consensus is growing that MgB2 
is describable within the conventional BCS-electron- 
phonon framework. However, we have proposed in ref. 
[Q that instead MgB2 should be described by the model 
of hole superconductivity, just as the high Tc cuprates 
|lO| . The common elements in the two classes of materi- 
als are that conduction is dominated by carriers in nearly 
filled bands, i.e. of hole-like character, and that the 
carriers that drive superconductivity propagate in con- 
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ducting substructures that are highly negatively charged 
(planes). Differences in the behavior of the two classes 
of materials arise within the model of hole superconduc- 
tivity from the fact that they are in different parame- 
ter regimes: the high Tc cuprates are in a substantially 
stronger coupling regime (as indicated by the shorter co- 
herence lengths), particularly for low hole doping [ p8| . 
Compelling aspects of the theory of hole superconductiv- 
ity are that it can describe superconductivity in a wide 
range of coupling regimes, and that it could be a univer- 
sal theory of superconductivity for all materials . 

In the weak coupling regime, the predictions of the 
model of hole superconductivity are similar to those of 
conventional BCS theory, and hence to the predictions 
of weak coupling electron-phonon BCS theory. Hence 
experimental evidence for BCS behavior in e.g. temper- 
ature dependence of the gap ||l3|] , or in tunneling charac- 
teristics jT^-^lSt should not be taken to favor the electron- 
phonon model over the model of hole superconductivity. 
The isotope effect, conventionally assumed to favor the 
electron-phonon model, is also expected within the model 
of hole superconductivity , and hence should also not 
be used to differentiate between both models. An Eliash- 
berg analysis of fine structure in tunneling characteris- 
tics above the gap energy, that traditionally has been as- 
sumed to be the strongest proof for the electron-phonon 
mechanism, has not yet been performed for this material. 

Here we have focused on two properties that show 
a clear difference in the electron-phonon and the hole 
model. One is the hole doping dependence of the crit- 
ical temperature, which the hole model predicts to be 
much stronger than the electron-phonon model. These 
experiments have not yet been performed, and once ex- 
perimental results become available it will be possible to 
ascertain which of both models is favored. Of course it 
is possible that even if experiments show that supercon- 
ductivity is rapidly suppressed with hole doping, as the 
hole model predicts, electron-phonon theory may also ac- 
count for it if a rapid decrease of electron-phonon matrix 
elements or of the relevant phonon frequencies with hole 
doping is postulated to occur, or if a rapid increase in the 
Coulomb pseudopotential fj,* with hole doping is postu- 
lated to occur. If so, electron-phonon theory (e.g. in it's 
EP2 version) and the hole theory will become increas- 
ingly indistinguishable. 

The other property that shows a clear difference in the 
electron-phonon models and the hole model is the pres- 
sure dependence of Tc for uniaxial pressure that modifies 
the intra-plane B — B distances: the hole model pre- 
dicts a strong increase in Tc, and the electron-phonon 
model a (weaker) decrease in Tc- Again, once experimen- 
tal results become available it will be possible to decide 
which model is favored. However, here again it is pos- 
sible that even if experiments show that Tc is strongly 
enhanced by reduction of the intra-plane B — B distance, 
electron-phonon theory could account for it if a concomit- 
tant increase of electron-phonon matrix elements or of 
the relevant phonon frequencies is postulated to occur or 



a concomittant decrease of Coulomb pseudopotential is 
postulated to occur. 

Assuming experimental results will show a rapid de- 
crease of Tc with hole doping as predicted by our the- 
ory, it is interesting to examine how the parameters in 
electron-phonon models would have to change to account 
for such a drop, given the band structure results for the 
density of states shown in Figure 1. Figure 11 shows three 
possible scenarios: (1) the behavior required for the aver- 
age of the square of the electron-phonon matrix element 
< > versus hole doping of the B atoms, assuming 
the phonon frequencies of the relevant phonons and the 
Coulomb pseudopotential ^* stay constant; (2) the be- 
havior required for the phonon frequencies of the relevant 
phonons < ojiog >, assuming the electron-phonon matrix 
elements and /j,* stay constant, and (3) the behavior re- 
quired for the Coulomb pseudopotential ^* , assuming the 
electron-phonon matrix elements and phonon frequencies 
of the relevant phonons stay constant. It can be seen that 
in all cases a rather rapid variation of parameters with 
hole doping is required. Of course a suitable combina- 
tion of decrease in electron-phonon matrix elements and 
phonon frequencies and increase in /j,* could also account 
for such behavior. It should be stressed that such rapid 
variations of electron-phonon matrix elements, phonon 
frequencies and/or /i* with hole doping have so far not 
been predicted by the electron-phonon models [^j^ . 

The range of hole doping where superconductivity oc- 
curs in our model is not strongly dependent on the pa- 
rameters in the model for a wide range of parameters, as 
was shown in Figures 3 and 4. Hence the prediction that 
superconductivity should only occur in a narrow range of 
doping around MgB2 is a strong prediction of the model. 
If the band structure results for the position of the Fermi 
level are correct, it implies that MgB2 is somewhat over- 
doped in our model, and hence doping with electrons 
should increase Tc. This is in apparent contradiction 
with experimental results for Af(7i_2;^^a;S2 18|. Possible 
explanations for the discrepancy may be problems with 
sample quality, or that the lattice constants change with 
increasing Al content. 

We have also examined here the predictions of the 
model of hole superconductivity for various observables 
in a range of parameters that appears to be appropri- 
ate for MgB2- Because this is a weak coupling regime, 
most properties are found to be very close to conven- 
tional BCS behavior. If the appropriate bandwidth is 
rather small the universal asymmetry in tunneling pre- 
dicted by the theory becomes of appreciable magnitude. 
We found the penetration depth decreasing monotoni- 
cally with hole doping (assuming the clean limit), which 
is the same qualitative behavior seen in high Tc cuprates 
p9[ . The coherence length was found to increase mono- 
tonically with hole doping, which is also seen in high Tc 
materials. If experiments confirm these predictions they 
will support the proposed commonality in the physics of 
superconductivity in the cuprates and in MgB2-deriyed 
compounds. Given this behavior, in the absence of dis- 
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order a crossover from type II to type I behavior with 
hole doping should eventually occur. While disorder is 
likely to prevent this in alloys, we expect that if similar 
stochiometric compounds with larger hole content than 
MgB2 are found they will have a smaller Ginzburg Lan- 
dau parameter k than MgB2 {n ~ 26) |25| and possibly 
even be type I. 

Superconducting properties as function of hole dop- 
ing have not yet been discussed within electron-phonon 
models. We note however that the strong increase of 
density of states expected with hole doping suggests that 
electron-phonon theory may describe a crossover to a 
stronger coupling regime with hole doping, i.e. decreas- 
ing coherence length and increasing penetration depth. 
We stress that this would be in qualitative disagreement 
with our predictions. 

In future work we will examine the predictions of our 
model taking into account the presence of two different 
hole bands, and the effect of the anisotropic band struc- 
ture, in order to calculate observables in different direc- 
tions which will be of interest once experimental results 
in single crystals become available. 
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FIG. 1. Comparison of the predictions for variation of the 
critical temperature with hole doping in the model of hole 
superconductivity (full line) and the electron-phonon model 
(dashed line). Here and in the following figures, is the 
average hole content per Boron atom; the total hole doping 
per unit cell is approximately three times larger. The results 
for electron-phonon theory were obtained assuming constant 
electron-phonon matrix elements and phonon frequencies, and 
using the density of states values obtained in the band struc- 
ture calculation in Ref. 2. 



FIG. 2. Comparison of the predictions for variation of the 
critical temperature with in-plane B-B distance in the model 
of hole superconductivity (full lines) and the electron-phonon 
model (dashed lines). The bandwidth and density of states 
are assumed to change by 5%. Again, the results for elec- 
tron-phonon theory were obtained assuming constant elec- 
tron-phonon matrix elements and phonon frequencies under 
compression. 



FIG. 3. Tc versus hole concentration for bandwidth 
D = 5eV and various values of Coulomb interaction parame- 
ters, given in the Figure in eV. The values for the correlated 
hopping parameter used for the four cases shown are in eV , 
K = 2.97,4.115,5.135,4.134 respectively, in the order given 
in the Figure label. 
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FIG. 4. Tc versus hole concentration for on-site Coulomb 
repulsion U = 5eV W = 0, and various values of the band- 
width D, given in the Figure in eV. The values for the cor- 
related hopping parameter used for the three cases shown are 
in eV, K = 2.97, 2.496, 1.963 respectively, in the order given 
in the Figure label. 

FIG. 5. Coherence length (in units of the lattice spacing) 
versus doping for the three sets of parameters of Figure 4. A 
lattice spacing corresponds to 2.2 A. 

FIG. 6. London penetration depth versus doping for the 
three sets of parameters of Figure 4. 

FIG. 7. Gap versus hole doping for the three sets of par- 
rameters given in Figure 4. 

FIG. 8. Specific heat jump versus doping for the three sets 
of parameters of Figure 4. 

FIG. 9. Temperature dependence of specific heat (a) and 
of energy gap (b) for the parameters of Figure 1 and doping 
nn = 0.065. 

FIG. 10. Tunneling characteristics for the parameters of 
Figure 4 with D = leV, doping Uh = 0.065, and tempera- 
tures T/Tc = 0.99,0.9,0.3,0.1. Note the higher peak when 
the sample is negatively biased. For larger bandwidths the 
magnitude of the asymmetry decreases. 

FIG. 11. Variation of the electron-phonon matrix el- 
ement < > (solid line), the average phonon frequency 
< uJiog > (dashed line), and the Coulomb pseudopotential 
H* (dash-dotted line) required, if Tc is given by the elec- 
tron-phonon model Eq. (7), to yield the rapid drop of Tc with 
hole doping predicted by the model of hole superconductivity. 
In each case it is assumed that the other parameters are fixed 
at their values for MgB2. Density of states values obtained 
from Ref. 2 are used. 
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